To the positive miscut influence 
on the crystal coUimation efficiency 

V. V. Tikhomirov, A. 1. SytoAfl 

Research Institute for Nuclear Problems, 
Belarus State University, Bobruiskaya 11, 220030 Minsk, Belarus 
(Dated: September 27, 2011) 

Abstract 

The paper concerns the crystal based coUimation suggested to upgrade the Large Hadron Cohider 
cohimation system. The issue of colhmation efficiency dependence on the muscut angle character- 
izing nonpar allelity of the channeling planes and crystal surface is mainly addressed. It is shown 
for the first time that even the preferable positive miscut could severely deteriorate the channeling 
collimation efficiency in the crystal colhmation UA9 experiment. We demonstrate that the positive 
miscut influence can increase the nuclear reaction rate in the perfectly aligned crystal collimator by 
a factor of 4.5. We also discuss the possible miscut influence on the future LHC crystal collimation 
system performance as well as suggest simple estimates for the beam diffusion step, average impact 
parameter of particle collisions with the collimator and angular divergence of the colliding particle 
beam portion. 
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I. INTRODUCTION 



Crystal based collimation was proposed to facilitate the beam halo cleaning at large 
accelerators long ago. Its application to the LHC upgrade becomes more and more topical 
The basic idea is to use a bent crystal in channeling mode to deflect halo particles 
w relatively large angles to high impact parameters of particle collisions with an absorber 

m 

If the first particle collision with the crystal collimator occurs at sufficiently small particle 
incidence angle w.r.t. the crystal planes (at pure alignment, 6c = 0), the probability of 
particle capture into the channeling regime reaches its maximum However even a small 
nonparallelity of the lateral crystal surface with atomic planes, characterized by the muscut 
angle 6^,, is able to severely disturb the motion of particles hitting the crystal with small 
impact parameters. In particular, if the miscut angle is negative, the channeling motion 
can be interrupted before the particle reaches the exit crystal face [sj. Since a considerable 
number of such particles will not reach an absorber fast, the negative miscut is recommended 
to be avoided |5|. By this reason the positive one (see Fig. 1) was chosen for the recent 
UA9 experiments 2j aimed to demonstrate the viability of crystal collimation. However the 
nuclear reaction rate in the perfectly aligned crystal collimator about five time exceeding the 
theoretically predicted Q value was observed. In this paper we for the first time investigate 





FIG. 1: A crystal with positive miscut angle before (left) and after (right) bending with angle 
if. 6c (positive), 6m (positive) and 6s = 6c — 6m (negative) are, respectively, the crystal plane 
misalignment angle, miscat angle and crystal surface misorientation angle, all measured in the 
direction of crystal bending, at that the angles 6c and 6s - from the z axis, parallel to the velocity 
of the particle just touching the crystal, and the angle 6m - from the crystal surface direction. 
Particles, moving from the left to the right with small impact parameters, enter the crystal through 
the lateral (upper) crystal surface. 
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the influence of positive miscut on collimation efficiency and demonstrate that it gives rise 
to up to 4.5 time increase of the nuclear reaction rate in the collimator. We also predict the 
low influence of the positive miscut on the efficiency of the future LHC crystal collimation 
system but first suggest simple estimates for the beam diffusion step, impact parameter 
of particle collisions with the collimator and angular divergence of the part of the beam 
particles colliding with the latter for the first time. 

II. PARTICLE DIFFUSION IN THE ACCELERATOR RING 

When positively charged particles strike a modestly bent crystal moving strictly parallel 
to its planes, they are captured into the regime of stable channeling motion with a proba- 
bility of 80-85% [6]. However if the incident particle beam possesses an angular divergence, 
the "channeling probability" decreases by the value APch oc (t?^), where (t?^) is the mean 
square of the incident beam divergence angle, assumed here to be considerably smaller than 
the critical channeling angle -^ch- This decrease remains negligible only if (^9^) < O.OW'^i^. 
Fortunately, the angular divergence of the beam portion striking a primary collimator first 
time often satisfies this rigorous condition. However if a particle has not been captured at 
the first passage through the collimator, its deflection can reach significant values. This 
actually pertains to the 15-20% of particles escaping channeling at the crystal penetration 
through its normal entrance (transverse to the beam) face. Particles entering a crystal 
with negative miscut at small enough impact parameters are angularly dispersed even 
stronger. Because of this and also since a miscut can not be avoided in practice, the positive 
miscut crystal orientation is commonly preferred [5|. In this case, however, the particles 
with the small enough impact parameters enter the crystal through its lateral face. Most 
of them avoid capture even in the case of pure crystal alignment. The uncaptured particles 
are scattered nearly the same way as the ones in amorphous matter, acquiring the average 
deflection angle squared proportional to the length Az of particle path through the crystal. 
Besides the angles of miscut 9^, crystal bending ip = l/R and crystal plane misalignment 
6c, the Az value depends on the particle impact parameter A with the crystal collimator, 
the value of which needs special consideration. 

Having limited access to the parameters of particle motion in the accelerators, not men- 
tioning the tools allowing to simulate a number of them, we need to suggest simple estimates 
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TABLE I: Beam diffusion parameter. 
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for the former. To start with, recall that particle collisions with crystal collimator originate 
from the particle diffusion in the radial beam direction caused by intra beam collisions, 
scattering by residual gas, elastic scattering at the interaction point, etc. {4]. Since a joint 
description of all these processes is hardly available, we will proceed from a simple estimate 
based on the accelerator beam lifetime r, particle revolution period T, r.m.s. beam radius 
a and collimator radial coordinate pc- We will also assume, as in j^], for simplicity, that the 
beam is axially symmetric and possesses normal distribution 



dN Np ( p2 



in particle radial number density integrated over the particle revolution period T. Here N is 
a total number of particles in the ring. Introducing a diffusion step 5, an average increase 
of radial coordinate acquired during one revolution by particles reaching the collimator 
position, one can express the particle loss rate in two ways: 

dN _N _ (dN\ 5 

coming directly to an estimate 

Table I illustrates Eq. (3) application to the cases of both the UA9 experiment and IR7 beta 
collimation region at the LHC [7]. Note that the exponential dependence on the collimator 
aperture squared leads to a really drastic 6 difference in the cases of collimation dedicated 
UA9 experiment with low intensity beam and pc — 3.5cr and the intensive LHC beam and 
Pr = 6cr. 
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III. PARTICLE IMPACT PARAMETER AND DEFLECTION ANGLE 



To simulate both the particle impact parameter and angular deflection at the moment 
of the first collision with the collimator we will proceed from the usual pseudoharmonic 
representation 

x{ip) = xqcosiP (4) 

of the betatron oscillations. Here ip and xq — yje^ are, respectively, the oscillation phase 
and amplitude, the latter of which is determined by the beam emittance e and accelerator 
beta function ^. Particle collisions with the collimator become possible since the amplitude 

Xq reaches the transverse collimator coordinate Xc = Pc (from here on we consider betatron 
oscillations in some transverse plane parallel to x axis). Collision really occurs if x{ip) > Xc 
or 1-01 < ip{xo), where 

iP{xq) = arccos ( — ) - J'^{xo - Xc)/xc (5) 

(we assumed that x(il^) — Xc ^ Xc). The particle direction of incidence on the crystal will 
be described by the angle 



simp —cosip 

2 as 



(6) 



of the velocity deviation from the direction of motion of the particle just touching the 
collimator having xq = Xc- 

The process of increase of the betatron oscillation amplitude was simulated using the 
formula 

Xo{n) ^Xo{n-l) + 2SCi, (7) 

where n = 1,2, .. is the number of particle revolution in the ring after the moment when 

a;o(0) = Xc- G (as well as .^2 and ^3 below) are random numbers uniformly distributed 
through the interval (0, 1). Note that in average {xQ{n)) = (xo(n — 1)) + S, in agreement 
with the S definition. 

In absence of phase correlations between betatron oscillations on different revolution 
periods the phase can be sampled by the formula ip — (2^2 — l)?!". A collision occurs with 
the probability 

p„ = ^ti^ c (8) 

TT Xr 
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at the n-th revolution if 



< ^(xo(n)). 



(9) 



If, otherwise, > ilj{xQ{n)), no colhsion occurs and one should continue the simulations 
with n = n + 1 and so on. Since the cumulative collision probability 



N 



n=N 
n=l 



Svr V Xr 



(10) 



increases faster and faster with the revolution number N, the "collision condition" (9) in- 
evitably becomes fulfilled at some revolution with some random values of xo{N) and 
allowing to evaluate both the collision coordinate (4) and angle (6) of particle deflection at 
the moment of collision. Particle distributions in the impact parameter A = x — Xc and 
deflection angle 6, simulated for the UA9 and LHC cases, are represented in Figs. 3-6. A 
three order in value difference in the impact parameter in the UA9 and LHC cases is directly 
related to that in the diffusion step - see Table I. 

Assuming that the colhsions, naturally, occur at P/v ~ 1, Eq. (10) can be reversed to 
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FIG. 2: A particle entering the crystal through the lateral surface. The crystals extends from 
z = to z = /, Xc is the collimator radial coordinate, {Xf> = Xs{ze),Ze) is the particle enter point 
while (xs(0),0) is the left upper corner of the crystal, 9 and Og are the deflection particle angle 
and crystal surface misalignment one, respectively, Az is the length of particle trajectory inside 
the crystal, 6 and A are, respectively, the particle diffusion step and impact parameter. 
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FIG. 3: Particle distribution in impact parameter for the UA9 case. 

estimate the typical revolution number before the collision 

^^(37rP^3/x7^ 3/^^ (11) 
2 \ S V (5 ' 

the average impact parameter 

(A) = Pn[xo{n)cos{^) - Xe] ^ ^ ^ l.^fj' (12) 

n=l 

and absolute value of the deflection angle 

m) = Pn' E PnXo{n)s^n{^)/f3 ^ ~ 0.8 (13) 

n=i 8 /:( 

The right hand sides of Eqs. 11-13 contain numerical factors found from the simulations 

' 1 ' 1 ' 1 ' 1 — 
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FIG. 4: The same for the LHC case. 
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FIG. 5: Particle distribution in deflection angle for the UA9 case. 



and determining some effective collision probabilities P/v, P^, ~ 1 which should be 
considered as the parameters compensating a slight model inconsistency. Despite the latter 
Eqs. 11-13 provide a sufficient base to compare the conditions of the UA9 crystal coUimation 



'ormance as well as to estimate the 
8| . Fig. 7 again demonstrates that 



experiment with the possible LHC crystal collimator per: 
perspectives of possible crystal coUimation development 
the radical difference in diffusion steps results in the drastic difference in average impact 
parameters in the UA9 and LHC cases. 
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FIG. 6: The same for the LHC case. 
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FIG. 7: Average impact parameter vs average beam diffusion step for tlie SPS UA9 (upper curve) 
and tlie LHC (lower one). Solid parts mark the actual parameter regions. 

IV. COLLIMATION EFFICIENCY FOR THE UA9 EXPERIMENTS 

A simulated value of the impact parameter A can be directly used to evaluate both the 
entrance transverse coordinate (4) and angle (6), the knowledge of both of which is necessary 
to simulate the particle trajectory inside the crystal in order to obtain the angle of particle 
deflection by the latter. Recall that we measure the angles from the direction of motion of 
the particle just touching the collimator at its maximum displacement x{iIj = 0) = Xq = Xc 
from the beam axis - see Fig. 2. Note de bene esse that the chosen zero direction forms the 
angle 

dx{ip = ^) 1 Fedp Xc , . 

w.r.t. the beam axis. Here s is the beam longitudinal coordinate and a = —dP/ds/2 is the 
conventional Twiss parameter. In the ideal case a crystal has no miscut and its planes form 
zero angle w.r.t. the chosen zero angle direction. Particles neither enter the crystal through 
its lateral surface no leave the one through it if they are channeled in this case. 

The real situation is complicated by the inevitable presence of both crystal miscut and 
crystal plane misalignment at the entrance surface, characterized by the angles dm and 6c, 
respectively. The crystal misalignment angle is assumed to be positive if the planes are 
rotated in the direction of crystal bending. If one determines the miscut angle as the one of 
crystal plane rotation in the direction of crystal bending w.r.t. the crystal lateral surface, 
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the misorientation angle of the latter, measured from the same zero angle direction, will be 
equal to = ^s(O) — 9c — Om- If the crystal is bent with radius R, the surface tangential 
direction will vary like 9s{z) = 9so + z/R with the longitudinal coordinate z s — Sc- A 
behavior of the surface coordinate 



also measured in the crystal bending direction, considerably differs if > and < 
and, in the latter case, if |6'so| > and l^^ol < V, where (p — l/R is the bending angle of 
the crystal with length /. Namely, if < a particle can enter the crystal through the 
lateral surface and, if —(p < 9so < 0, also leave it through the latter. On the opposite, if 
^so > particles always enter the crystal through the entrance face, while leave it either 
through the lateral or exit ones. In all the cases the actual situation is determined by the 
impact parameter A, the random nature of which allows for diverse trajectory types at any 
choice of 9s and (p. To get simple formulae for all possible situations we first determined 
the minimal crystal surface coordinate Xmin = Xc- After the Monte Carlo sampling of the 
impact parameter value A = x{iIj) — Xc corresponding transverse collision coordinate x{iIj) 
was used to evaluate the longitudinal one Ze from the equation Xs{ze) — x{%Ij). Then both 
the particle entrance point coordinates {x{'il)),Ze) and initial deflection angle (6) were used 
as the initial conditions for its trajectory simulation inside the crystal. A possibility to leave 
the crystal through the lateral surface at some z < I was permanently monitored using Eq. 
(15). The particle transverse coordinate and deflection angle at the exit together with the 
crystal position in the ring became the initial conditions for the particle motion simulation 
in the accelerator ring, for which the simplest model of betatron motion was applied. 

However since the simulation of thousands of trajectories for tens of different miscut angle 
and diffusion step values takes considerable time, more rationally was flrst to use a simplifled 
fast approach allowing to elaborate a general view on the influence of positive miscut on 
the coUimation efficiency. The idea originates from the mentioned proportionality of the 
decrease of the channeling efficiency at the second crystal penetration to the squared angle 
of multiple scattering of particles entering the crystal first through the lateral surface. Since 
the latter, in turn, is proportional to the scattering length A2;, one can conclude that simply 



Xs{z) = Xs{0) + f 9s{z)dz = Xs{0) + 9soZ + z'^/2R, 
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FIG. 8: Measure in centimeters average length (Az) of scattering of particles entering the crystal 
through the lateral crystal surface vs both miscut angle and diffusion step at perfect crystal plane 
alignment. 



and reduce the issue to the analysis of the behavior of the averaged length (Az) of the first 
pass through the crystal of the particles entering the latter exclusively crossing its lateral 
surface. Fig. 8 illustrates the simulated behavior of (Az) in the typical UA9 case of / = 2mm 
and if = 1/ R = 150 firad. Surprisingly, the simulations unambiguously point to the region 
9m ~ lOO^rad and 5 ~ 1 A of the UA9 experiment parameters as to the one of the greatest 
possible miscut influence on the coUimation. 

Since the scattering length determines the decrease of the channeling probability and since 
nonchanneled particles induce more nuclear reactions that the channeled ones, Fig. 8 has to 
simultaneously reflect the behavior of the rate of the nuclear reactions induced in the crystal 
collimator. To illuminate the possible role of positive miscut in the UA9 experiment we, 
according to Ref. and Table I, had put 6 = 1 A and conducted more detail Monte Carlo 
simulations of the nuclear reactions in the miscut angle interval —300 fir ad < 6m ^ 300 fir ad 
taking now into detail consideration also the particle motion in the crystal collimator. The 
dependence obtained (see Fig. 9) demonstrates an evident agreement with that of (Az) 
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FIG. 9: Probability of nuclear reactions in the crystal collimator vs miscut angle at perfect crystal 
plane alignment. 

along the vertical (red) line drawn at 5 = 1 A in Fig. 8, confirming thus the strong influence 
of the positive miscut on the collimation process. In principle, positive miscut causes even 
slightly larger increase in nuclear reaction rate (the right peak) than the negative one (the 
left one). At this the increase of the reaction probability caused by the positive miscut with 
dm — 125firad reaches 8.6/1.9 ^ 4.5. Thus, the miscut influence, for sure, should be taken 
into consideration for the full interpretation of the UA9 experiments j^. 



V. MISCUT INFLUENCE AT THE LHC 



For the future possible application at the LHC it is important is to clarify how the 
deteriorating miscut influence can be avoided. A joint consideration of the particle motion 
in both the ring and the crystal results in the encouraging conclusion that the observed 
undesirable increase in nuclear reaction rate can be easily avoided in both UA9 and LHC 
cases. In fact, some of the conditions of the UA9 experiment prove to be practically optimal 
for the demonstration of the maximum miscut role. The point is a perfect matching of 
the average impact parameter (12) (A) ~ 0.039/im with the width Xs{0) — Xc = 6'^R/2 ~ 
0.067/im of the impact parameter region allowing particle entrance through the lateral crystal 
surface. This matching made possible both the lateral entrance of the majority of particles 
and their relatively continuous path inside the crystal, the average value of which (A^;) ^ 
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1.2mm exceeds a half of the crystal length / = 2mm. It is namely the nearly "amorphous" 
scattering at such a length which gave the origin to the angular dispersion of particle beam 
causing the decrease of the capture probability to the channeling regime at their subsequent 
passages through the crystal collimator. 

At least two ways to decrease the miscut role by fulfilling the condition A ^ Xs{0) —Xc< 
0ml can be readily suggested. The most evident, though, probably, more difficult, is to lessen 
the miscut angle down to about lOyurac?, as Figs. 8 and 9 suggest. The second is to increase 
the collision parameter (12) (A) oc 5"^^^ by means of beam diffusion acceleration. While 
the diffusion step could be nearly freely chosen in the coUimation UA9 experiment, the 
actual set of the LHC parameters solves this problem automatically. Indeed, if i? = 100m 
and / = 4mm, one obtains Xs(0) — — 0.32/im, or more than a hundred times less than 
(A) ^ 43yum ~ 130(a;s(0) — Xc) without special measures. Thus, only a negligible portion of 
the LHC protons will enter the crystal collimator through the lateral crystal surface even at 
the typical miscut angles of Om ~ lOO/irarf. 

It also should be noted that despite the relatively large value of the diffusion step 5, the 
angular divergence (13) of the colliding beam portion, as Fig. 6 demonstrates, is low enough 
to provide the probability of capture into the regime of channeling motion comparable to 
the maximum one. Nevertheless some decrease in divergence remains desirable. The sharp 
dependence (3) of the diffusion step on the collimator aperture Xc/u allows to decrease 6 
by means of a slight decrease of the latter. At this, if the divergence of the colliding beam 
portion is decreased by several times, it will become possible to sharply rise the probability 
of particle capture into the channeling regime up to 99% by the method of the crystal cut 



In conclusion, the positive miscut influence indeed could increase the nuclear reaction 
probability in the crystal collimator up to 4.5 times. Nevertheless if the crystal collimator 
system based on the channeling particle deflection is realized at the LHC, its functioning 
will not be considerably disturbed by the influence of crystal miscut. In addition, the 
performance of the crystal collimator can be drastically improved by the method ^ of the 
crystal cut. 
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